Energy landscape and shear modulus of interlayer Josephson vortex systems 



o 
o 

(N 



(N 



Yoshihiko NonomurgQ and Xiao Hu 

Computational Materials Science Center, National Institute for Materials Science, Tsukuba, Ibaraki 305-0047, Japan 

(Dated: February 6, 2008) 

The ground state of interlayer Josephson vortex systems is investigated on the basis of a simpli- 
fied Lawrence-Doniach model in which spatial dependence of the gauge field and the amplitude of 
superconducting order parameter is not taken into account. Energy landscape is drawn with re- 
spect to the in-plane field, the period of insulating layers including Josephson vortices, and the shift 
from the aligned vortex lattice. The energy landscape has a multi- valley structure and ground-state 
configurations correspond to bifurcation points of the valleys. In the high-field region, the shear 
modulus becomes independent of field and its anisotropy dependence is given by C66 oc 7~ 4 . 

PACS numbers: 74.25.Qt 
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Introduction. Although interlayer Josephson vortex 
systems in cuprate high-T c superconductors have been in- 
tensively studied, their phase diagrams have not been es- 
tablished yet even in the ground state. For high in-plane 
fields, Josephson vortices penetrate into every insulating 
layer, and the ground state is given by an elongated tri- 
angular lattice aligned along the superconducting layers. 
As the field decreases, such an aligned vortex lattice be- 
comes unstable owing to the shear instability , and the 
shearing angle of vortex lattices increases continuously 
0- As the field further decreases, the rotated vortex lat- 
tices have been considered to be the ground state 0,0,0 • 
The studies mentioned above (except for Ref. Q) were 
based on the London model, and effects of the layered 
structure were included only as geometrical constraint. 

This layered structure can be directly taken into ac- 
count by using the Lawrence-Doniach model. Bulaevskii 
and Clem Q first introduced this model in order to in- 
vestigate interlayer Josephson vortex systems for high 
fields. Ichioka [?J calculated the energies of the aligned 
vortex lattices including vacant insulating layers. Sta- 
bility of some rotated-vortex-lattice configurations was 
first pointed out by Hu and Tachiki in the frustrated 
XY model, and Ikeda [9| reported similar results in the 
Lawrence-Doniach model within the lowest-Landau-level 
approximation. Quite recently, Koshelev |2j] systemati- 
cally evaluated the ground-state phase diagram (includ- 
ing the sheared and rotated vortex lattices) based on a 
simplified Lawrence-Doniach model, in which spatial de- 
pendence of the gauge field and the amplitude of super- 
conducting order parameter is neglected. 

The ground state of this simplified model is completely 
given by the in-plane field h, the period N between insu- 
lating layers including Josephson vortices, and the shift 
A of the vortex lattice from the aligned one. In Ref. 
, A is fixed at fractional numbers which correspond to 
the rigid rotated vortex lattices. In the present study, we 
draw the full energy landscape for the first time, and find 
that the energy landscape has a multi-valley structure, 



and that the ground state for low fields shows systematic 
deviation from the rigid rotated vortex lattices. We also 
find that field dependence of the shear modulus is quite 
different from previous estimation based on the London 
theory pj in the high-field region. 

Formulation. When spatial dependence of the gauge 
field and the amplitude of superconducting order param- 
eter is neglected, the Lawrence-Doniach free energy func- 
tional is expressed by the phase component <j) n [3] as 
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where <po is the flux quantum, A a fc is the penetration 
depth in the ab plane, 7 is the anisotropy parameter, 
and d is the interlayer distance. The distance a between 
Josephson vortices in the same layer and the period N 
(see Fig. Ufa)) are related with the in-plane field along 
the x axis B x as cf>o = adNB x . Since the ground state 
should be uniform along the field direction, the free en- 
ergy functional per unit volume is expressed as 
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where the following normalized quantities are introduced: 
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By minimizing the free energy functional © with respect 
to cf> n , we have 
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FIG. 1: Definition of (a) the unit cell of a vortex lattice 
and (b) the labeling of a rigid rotated vortex lattice (here 
(n,m,N) — (2, 1, 1)). Scaled coordinate y = y/^d) is used. 



FIG. 2: (a) Energy function —g versus the shift from the 
aligned vortex lattice A and the normalized in-plane field h. 
Each color stands for the energy function for each TV [color on- 
line] . The dark [blue] solid lines denote the ground state, and 
the pale [green] solid lines stand for the lowest-energy state 
(but not the ground state) for each N. The rotated-lattice in- 
dex (n, m, N) for each ground state configuration is displayed, 
(b) Top view of the above figure [gif figure appended] . 
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This differential equation is solved in the unit cell a x Nd 
(see Fig.^a)) with appropriate quasi periodic boundary 
conditions [1C|. In order to search the ground state, the 
shift from the aligned vortex lattice A defined in Fig. 
^a) should also be swept as the in-plane field h. 

Another description of the vortex lattice in layered 
material is shown in Fig. ^b), where the vector a be- 
tween the nearest-neighbor vortices in the same layer 
and the unit vectors of the triangular vortex lattice ei 
and e2 are combined with a pair of integers n and to as 
a = nei + me^. This labeling is based on an equilateral 
triangular vortex lattice, which we call the "rigid rotated 
vortex lattice" in the present article. For each (n,m), 
there is a corresponding fractional value of A. However, 
it is not trivial that any rigid rotated vortex lattice can 
be a ground state. As will be revealed below, it is not the 
case. Hereafter the energy function o = log h— 1.432 
is used instead of u, following Ref. Q|. 

Results. The energy landscape of eq. © is shown in 
Fig. EJa), where —g is taken instead of g because the 
valley structure is hard to visualize. The ridges in this 
figure correspond to the valleys in the actual energy land- 
scape. Energy functions for various iV's are drawn in the 
same figure, and different colors correspond to different 
TV's. The ground state denoted by dark solid lines of- 
ten jumps between two different values of A as drawn in 
dashed lines. The pale solid lines represent the lowest- 
energy state (but not the ground state) for each N. 

For high fields, the aligned (1, 1, 1) vortex lattice cor- 
responds to the ground state with A = 0. As the field 
decreases to h w 1.33, a sheared lattice with finite A 
becomes the ground state. The shearing angle increases 
continuously as the field decreases, and such continuous 
field dependence cannot be represented by the rotated- 
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FIG. 3: Ground-state configurations for various in-plane fields 
h with the corresponding rotated-lattice indices (n,m,N). 



lattice index. Although the ground state changes to the 
(1,0,2) aligned lattice at h ~ 0.99, the sheared vor- 
tex lattice state is still the lowest energy state in the 
N = 1 subspace, and it becomes the ground state again 
at h rs 0.80. The field dependence of A is quite small 
around h w 0.7, which corresponds to the simplest rigid 
rotated vortex lattice labeled with (2,1,1) and charac- 
terized by A r igid = 1/7. Since each ground-state con- 
figuration is separated in the energy landscape, these 
configurations can be labeled by the rotated-lattice in- 
dices (n, to, N), even though the parameter A may devi- 
ate from the fractional values A r ; g id corresponding to the 
rigid rotated vortex lattices. Typical ground-state config- 
urations are shown in Fig. [31 where all the rotated-lattice 
indices (n, m, N) for h > 0.15 are listed. 

Next, the energy landscape in the N = 1 subspace is 
shown in Fig.0] Here the dark or pale solid lines are de- 
fined similarly to the ones in Fig.Et a ), an d the paler solid 
lines denote the local energy minima (not the lowest- 
energy state in this subspace). Apparently, ground-state 
configurations of the rotated lattices are connected with 
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FIG. 4: Energy function — g versus A and h in the N — 1 
subspace [color online]. The dark [blue] or pale [green] solid 
lines are the same as in Fig. [51 a), and the paler [yellow] solid 
lines stand for the local energy minima (not the lowest-energy 
states), and its top view is displayed in the inset. The rotated- 
lattice indices are given for all the ground-state configurations 
in this subspace [gif figure appended]. 
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FIG. 5: Deviation of A from that of the rigid rotated vor- 
tex lattice (A r igid) versus normalized field h in the N = 1 
subspace, and it takes mi 111 mum fit h — hmin 

(dashed lines). 

The fields corresponding to the rigid rotated vortex lattices 
(frigid) are marked by arrows. In the inset, the deviation 
of A (| A r igid — A|) at h — /i m in is plotted versus that of h 
(/irigid — hixmx) together with the rotated-lattice indices. 



each other by lines of local energy minima, and such 
ground-state configurations correspond to the bifurcation 
points of local minima. Similar bifurcating structure is 
expected for N > 2 at much low fields. 

From precise analyses of A, we find that the rigid 
rotated vortex lattices expected from the London the- 
ory can never be the ground state. The deviation of A 
from the values corresponding to the rigid rotated lat- 
tices A rig id is plotted versus h in Fig. The inequal- 
ity A rig id — A > holds for n > to/2, and vice versa. 
The field corresponding to the rigid rotated vortex lat- 
tice ft. r igid is marked by arrows in the same figure, and 
|A r i g id — A| takes minima at h = /i m in (dashed lines). The 
deviation decreases as h decreases, and seems to vanish 
in the h — ► limit. In the inset, the deviation of A at 
h = h mm is plotted versus that of the field, /i r i g id — h mm . 

Finally, we observe the shear modulus cqq defined by 
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where <fi = adNB x and eq. (0J are substituted and pa- 
rameter dependence on h, N and A is only included in 
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FIG. 6: Field dependence of the scaled shear modulus C66 
defined in eq. ®. (a) For h>l, C66 converges to a constant 
value C w 39.4 for high fields, and the vanishing point of cee 
divides the aligned and sheared vortex lattices, (b) For h < 1, 
corresponding rotated-lattice indices are displayed. 



the scaled shear modulus Cqq. Field dependence of Cqq is 
displayed in Fig. El and this quantity converges to a con- 
stant value C w 39.4 in the h — ► oo limit. Namely, the 
shear modulus becomes field independent for high fields, 
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This behavior is in sharp contrast to the exponen- 
tially vanishing cqq of the London model and is 
also different from that of the elastic theory, cqq = 
[(j>oB x / (87rA a (,) 2 ]7 -3 [llj- Introducing a "saturated field" 
B defined below, eq. Q can be expressed similarly to the 
elastic theory as 
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and a = 4>o/(dB x ) = (2n 3 /C)jd « 1.57jd corresponds 
to the effective core size of a Josephson vortex. As the 
field decreases, cqq vanishes at the shear instability field 
of the aligned vortex lattice , increases again together 
with the increase of the shearing angle, and jumps when 
the rotated-lattice index changes. 

Discussion. In the London model, only the rigid vor- 
tex lattices can be the ground state by definition. There- 
fore, the sheared vortex lattices are specific to the mod- 
els including the layered structure explicitly such as the 
present model. We find that the rotated vortex lattices 
with Arigid cannot be the ground state even at the fields 
corresponding to the rigid rotated lattices, which is an- 
other nontrivial effects of the layered structure. 

It would be interesting to investigate how the above 
description is altered by various kinds of fluctuations. 
Hu et al. jl2j] studied the phase diagram of interlayer 
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Josephson vortex systems in the vicinity of the melting 
temperature using the density functional theory. They 
found the vortex smectic state around the (n, m, N) = 
(1, 0, 2) state, while no sheared vortex lattice states were 
observed between the (1,1,1) and (1,0,2) states. This 
finding may suggest that the sheared vortex lattice state 
is affected very much by thermal fluctuations. Study on 
temperature dependence of the deviation from the rigid 
rotated lattices would be an interesting future problem. 

The screening effect is also studied in the present 
framework, and we find that the rotated vortex lattices 
are almost unchanged, while the sheared vortex lattices 
or the N > 2 aligned vortex lattices are affected very 
much. These facts suggest that the rotated vortex lattice 
structure is robust against the screenin g eff ect. Details 
of this study will be reported elsewhere jl3| . 

Summary. Ground state of interlayer Josephson vor- 
tex systems is investigated on the basis of the simplified 
Lawrence-Doniach model, in which spatial dependence of 
the gauge field and the amplitude of superconducting or- 
der parameter is neglected. The free energy functional is 
minimized in the unit cell with one Josephson vortex. We 
draw the energy landscape with respect to the in-plane 



field, the period of insulating layers including Josephson 
vortices, and the shift from the aligned lattice. 

As the in-plane field decreases, the aligned vortex lat- 
tice along superconducting layers changes to the sheared 
vortex lattice and then to the rotated vortex lattices, 
in which the shift from the aligned lattice is approxi- 
mately given by that of the corresponding rigid rotated 
lattices. Systematic deviation from the rigid configura- 
tions is clarified for the first time. Ground-state config- 
urations of the rotated lattices are connected with each 
other by local minima in the energy landscape to form a 
multi-valley structure, and such ground-state configura- 
tions correspond to the bifurcation points of the valleys. 

In the high-field region, the shear modulus becomes in- 
dependent of field and its anisotropy dependence is given 
by C66 oc 7~ 4 , quite different from the London theory in 
which C66 decays exponentially with field. This quantity 
vanishes at the onset of the sheared vortex lattice, and 
jumps when the rotated-lattice index changes. 
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